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ABSTRACT 

wo: 

The non-leptonic two body decays of D mesons are studied in the framework 
of an improved factorization approximation. The final state interaction ef- 
\Q ! fects are taken into account assuming them dominated by nearby resonances. 

The agreement with experimental data is quite good. CP-violating asym- 
metries are predicted. 

\o ' 

^ ! 1. Introduction 

00 ■ 

A theoretical description of exclusive non-leptonic decays of charmed mesons based 
on first principles has not yet been achieved. Although the short- distance effects due to 
hard gluon exchange can be resummed and the effective hamiitonian has been constructed 
at next-to-leading order ,0 the evaluation of its matrix elements requires non-perturbative 
techniques. In this respect, a classical analysis based on QCD sum rules has been presented 
in three papers by Blok and Shifman.H However only the general trends were reproduced 
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by their analysis, while no agreement with current experimental data was obtained. More 
recently, Martinelli and collaboratorsi have proposed a procedure to study two-body non- 
leptonic weak decays in numerical simulation of lattice QCD. Since no numerical result 
o3 ■ have been obtained as yet, one has however to resort to models. 

We here present one such modelS based on the factorization approximation with an- 
nihilation terms and rescattering effects due to the resonances coupled to the final states, 
that has been rather successful to account for the experimental data about two-body de- 
cays of charged and neutral D mesons in PP and PV final states. This feature has made 
it possible for us to obtain reliable predictions for the related CP-violating asymmetries. 

2. Weak Decay Amplitudes 

The effective weak hamiitonian for Cabibbo allowed non-leptonic decays of charmed 
particles is given by 

H%= AS = % V ud V; s \C 2 s a 7m (1 - 75) c a v? ^ (1 - 75) dp + 



V2 



d u a 7m (1 - 75) c a s p 7^(1 -Ts)^ +h.c. . (1) 
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For AC = —AS processes the hamiltonian is obtained from the previous equation with 
the substitution s <-> d. The effective weak hamiltonian for Cabibbo-first-forbidden (CFF) 
non-leptonic decays reads 



H. 



AC=±1,AS=0 
eff 



G f 



-§{ v ud v; d [d + c 2 Qi] + Ks [d qi + c 2 q s 2 



- v ub v; b E CiQi} + h.c. 

i=3 

In Eq. (H) the operators are defined as 
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The operator Q\ in Eq. (|]) is obtained from (Qf) with the substitution (d — > s). 
a and (3 are colour indices (that we will omit in Eq. (|j)) and in the "penguin" operators 
q (q) is to be summed over all active flavors (u, d, s). 

For the Wilson coefficients we u sed the anomalous dimension matrices calculated at 
next-to-leading order.S Assuming A% IS = 300 MeV, at the scale /1 = 1.5 GeV for the 
"scheme independent prescription" (cfr Buras et al. in Ref. 0) we obtain G\ = —0.628, 
C 2 = 1.347, C 3 = 0.027, C 4 = -0.057, C 5 = 0.015, C 6 = -0.070. 

In the factorization approximation the matrix elements of H e ff are written in terms 
of matrix elements of currents, (V^) M = q' q and {A^,)^ = q' '7^75^. To evaluate these 
matrix elements,* we adopt the usual definition of the decay constants and the form 
factors.i The q 2 dependence of the involved form factors (fi(q 2 ), fo(q 2 ) and A (q 2 ) for 
PP and PV final states) is assumed to be dominated by the nearest resonances. The 
values /i(0) and /o(0) are fixed by using SU(3) symmetry and the semileptonic decay 
rate D° — > K~e + u. Since the data on D meson decays show large SU(3) breaking effects, 
in our fit we allowed a cs (= Aq~* s (0)) to be different by a cd = a cu (= A%~~* d ^ u \o)) and the 
values obtained by the fit are a cs = 0.59, and a cd = l.t 

In the W-exchange and annihilation terms, however, the large and time-like q 2 values 
needed, together with the suggested existence of resonances with masses near to the D- 
meson mass, make a prediction based on the lightest mass singularity unjustified. These 
terms depend on the matrix elements of current divergences between the vacuum and 
two-meson states. We write them, with the help of the equations of motion, in the way 



*The matrix elements of axial vector current and axial density between the vacuum and rj(r]') state are 
evaluated following the authors of Ref. 

direct QCD sum rule calculation^ of Ao(q 2 ) shows the q 2 dependence compatible with the pole, but 
a different SU(3)-breaking effects: a cs /a cu — 1.10 ± 0.05 at q 2 = 0. 
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indicated in the following examples: 



M 2 

< K-n + \d*(y e %\0 > = i(m s - m d ) < K-^ + \ sd |0 >= i (m s - m d ) — ^ W PP , 

JD 

< K-p + \d^(A d s ) tl \0>= t(m s + m d ) <R-p + \ s l5 d\0>= -(m s + m d )^ e*-p K W PV .(4) 

JD 

We use SU(3) symmetry for the matrix elements of scalar and pseudoscalar densities, and 
express all of them in terms of Wpp, Wpy, which are, in our approach, free parameters of 
the fit. Their magnitude turns out to be considerably larger than what one would obtain 
assuming form factors dominated by the pole of the lightest scalar or pseudoscalar meson, 
i.e. K*(U30) or K(A97). 

3. Final State Interaction Effects 

As far as final state interactions (FSI) are concerned, we assume that they are domi- 
nated by resonant contributions, and we neg lect the phase-shifts in exotic channels! In the 
mass region of pseudoscalar charmed particles there is evidence! for a J p = 0~ ix"(1830) 
(with T = 250 MeV and an observed decay to Kip) and a J p = 0~ 7f(1770) with T = 310 
MeV. These resonances have the right quantum numbers to construct an 0~ octect which 
can couple to PV final states. The hamiltonian is determined from charge conjugation 
and SU(3) symmetry. Analogously, the FSI for D — > PP decays should be dominated by 
the J p = + octect. There is evidence for the existence of K*° (with mass 1945±10±20 
MeV, width 201±34±79 MeV and 52±14% branching ratio in Kir). Unfortunately, no a 
isovector resonance has been observed up to now in the interesting mass region. However, 
we assume his existence and fix the mass with an equispacing formula. 

The description of rescattering effects for Cabibbo forbidden D° decays is complicated 
by the presence of a coupling with a yet unobserved / and f isoscalar resonances, which 
should be singlet-octet mixtures."'" In order to reduce the number of free parameters we 
assume the scalar resonances behave as the tensor mesons (J p = 2 + ), ^2(1270) and 
/2(1525) } This procedure relates the coupling constant in the strong hamiltonian with 
the mixing angle <fi between the singlet-octect part of fo and /q." 

In our model the FSI effect modifies the amplitudes, A w , in the following way: " 

A(D -> V h P k ) = A W (D -> V h P k ) + c hk [ex P (i5 8 ) - 1] J2 c h>k>M D -> V^P*,) (5) 

hlkl 

In Eq. (H) Chk are the normalized {J2 c hk = 1) couplings PPV and 

T(P) 

sin<5 8 exp(i<5 8 ) = — (6) 

where P is the resonance appropriate to the decay channel considered (n or K). 

*Note that for the PV final state the C symmetry forbids the coupling of the final state to a singlet part 

of an hypotethical ^-resonance. Thus, we need to fit only a single phase, <5. 

§ The f 2 is very weakly coupled to nir, and the fi has in turn a small coupling to KK. 

^For further details see Ref. □. 

"An analogous expression holds for PP final state. 
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4. Comparison with Experimental Data on Branching Ratios and Predicted 
CP-Violating Asymmetries 

Using our model to evaluate the weak amplitudes A w and modifying them with FSI 
effects we are able to write the rates for all Cabibbo-allowed two-body decays and for 
Cabibbo-forbidden D + and D£ as functions of eleven free parameters. Their values are 
fixed with a fit to the experimental data. The total Xt = 90-0 f° r 45 data (x 2 /dof = 2.6): 
25 are data points for Cabibbo-allowed decays (% 2 = 61.8), 12 for CFF D + and D+ 
(X 2 = 18.8), 4 for CFF D° -> PP { X 2 = 1.7), and 4 for CFF D° -> PV (V 2 = 7.7). In all 
cases the agreement with the data is quite good, but the Table 1 shows that our model 
needs an improvement in describing the decays in PV final state.** 

Table 1: Partial % 2 for each class of D decays. 



Decays 


# data 


x 2 


D° -> PP 


8 


8.44 


D+ -> PP 


5 


9.56 


Df -> PP 


4 


8.79 


D° -> PV 


12 


18.35 


D+ -> PV 


8 


29.55 


P>+ - PV 


8 


15.35 



It is well known that CP- violating effects show up in a decay process only if the decay 
amplitude is the sum of two different parts, whose phases are made of a weak (CKM) 
and a strong (final state interaction) contribution. If A\ and A2 denote the generic two 
weak amplitude contributing to the D — > f amplitude, the CP-violating asymmetry in 
the decay rates will be: 

2 gQMg) sinfe - S 1 ) 
a ° P |A 1 | 2 + |A 2 | 2 + 2K(A 1 ^) cos(5 2 -5i) U 

where 5, are strong phases. Now in the Cabibbo-first-forbidden D decays the penguin 
operators in Eq. (0) provide the different phases of the weak amplitudes A\ and A^. 

Having obtained a quite good description of the rates we may give reliable prediction 
on CP-violating asymmetries for D + and D° Cabibbo-first-forbidden decays. In Table 2 
we report CP-violating asymmetries in 10~ 3 unit; the central values are obtained choosing 
p = 0.2, r) = 0.3 and Vb — 0.040. The errors result from the variation of p and rj in the 
one-sigma region obtained in Ref. 0. 

As we can see in Table 2 large asymmetries (~ 10 -3 ) are predicted in our model; at 
this end large final state phase shifts and penguin contributions played the fundamental 
role. 

**This is the starting point for the authors of Ref. 0. 
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Table 2: CP- violating decay asymmetries for some D + and D Cabibbo forbidden decays. 



decay channel 


10 3 x a CP 


decay channel 


10 3 x a C p 


D+ _^ p o 7r + 


-1 17 ± 68 


D + -> A U A+ 


-0 51 ± 30 


D+ -> p+7T U 


+1 28 ± 74 




-1 43 ± 83 

_1_ • J- v/ — 1 — V 7 • v. J - 7 


no _v /r* "^ 

±7 — > i\ J\ 


n p.7 -\- n 'iQ 
— u.o ( it u.oy 


n° _i -7r°r)' 


-ru.yo x u.o/ 


L> — > K K 


— (J. 57 ± (J. 39 


7~l0 

D — > 7777 


n m i n on 

—0.50 ± 0.29 




+0.038 ± 0.022 


Z) — > T7??' 

/ / 


-0.28 ± 0.16 


D° -> A*~A+ 


+0.16 ±0.09 


D° -> 7T 7T 


+0.54 + 0.31 


£>° -> P+7T" 


+0.37 + 0.22 




-0.02 + 0.01 




-0.36 + 0.21 


D° -> A+A~ 


-0.13 + 0.08 






730 ^ 


+0.28 + 0.16 
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